Damped Bloch oscillations of cold atoms in optical lattices 
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The paper studies Bloch oscillations of cold neutral atoms in the optical lattice. The effect of 
spontaneous emission on the dynamics of the system is analyzed both analytically and numerically. 
The spontaneous emission is shown to cause (i) the decay of Bloch oscillations with the decrement 
given by the rate of spontaneous emission and (ii) the diffusive spreading of the atoms with a 
diffusion coefficient depending on both the rate of spontaneous emission and the Bloch frequency. 



I. INTRODUCTION 



In 1928 Bloch predicted that the coherent motion of crystal electrons in a static electric field should be oscillatory 
rather than uniform Q. Nevertheless, this phenomenon (known nowadays as Bloch oscillations), has never been 
observed in bulk crystals. Because of relaxation processes (scattering on lattice defects, phonons, etc.) the coherence 
of the system is destroyed before electrons complete one Bloch cycle. This obstacle has forced the researchers to look 
for other systems, where the Bloch period (which is inversely proportional to the magnitude of the static force and the 
lattice period) can be smaller than the characteristic relaxation time. This is realized in semiconductor superlattices, 
and in 1992 a direct observation of the current oscillations in semiconductor superlattices was reported 0, f|. It 
should be stressed, however, that in semiconductor superlattices the relaxation time only slightly exceeds the Bloch 
period and, thus, in practice one always meets the regime of damped Bloch oscillations. A detailed analysis of the 
relaxation processes in the semiconductor system appears to be a rather complicated problem and phenomcnological 
approaches are usually used to describe the decay of current oscillations. 
, Recently Bloch oscillations were observed in the system of cold atoms in an (accelerated) optical lattice 0. This 
system of atomic optics mimics the solid state system where the neutral atoms and standing laser wave play the roles 
""^ | of electrons and crystal lattice, respectively. Unlike the solid state systems, it may show both the regimes of undamped 
and damped oscillations. Indeed, the main relaxation process in the system "atom in a laser field" is the spontaneous 
emission of photons by the excited atom and it can be well controlled by choosing an appropriate detuning from the 
(— i . atomic resonance. Up to now only the regime of undamped oscillations (large detuning) has attracted the attention of 
researches [|[ |[ |[ @, |§| ||, [To| - However, keeping in mind an analogy with the semiconductor superlattices, the latter 
regime is also of considerable interest. In this present paper we study the process of decay of atomic Bloch oscillations 
due to the effect of spontaneous emission. In particular we address the question of the actual decay process and the 
, dependence of the decay rate on the system parameters. 

• i-H , 



II. THE MODEL AND APPROACH 

We shall take into account the spontaneous emission by using the standard approach, where its effect is characterized 
by a single constant 7 defined as the inverse radiative lifetime of the upper state. Namely, it is assumed that the 
diagonal and off-diagonal elements of the 2x2 density matrix of a fixed atom relax to their equilibrium values with 
the rates 7 and 7/2, respectively. When atom is in free space, its density matrix depends additionally on the position 
z of the atom and, in general case, the resulting master equation has a rather complicated form |ll], |l2|]. This equation 
can be considerably simplified in the limit of large detuning 6 3> Q (here f2 and 6 are Rabi frequency and detuning 
from the atomic resonance). In this case we can eliminate the internal structure of the atom and the master equation 
for the position density matrix has the form 



(i) 



(see, Ref. [13, 14l, for example) where 



H=^-^fcos 2 (k L z) + Fz (2) 
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is the effective Hamiltonian of the system and 

L u = cos(fciz) exp(mfciz) , \u\ < 1 (3) 

is the projection of the recoil operator on z-axis. The distribution P(u) of a random variable u in Eq. (0) is defined 
by the angle distribution for the momentum ftk.L = ft/c^n of the spontaneously emitted photons, which in the case of 
linearly polarized light is given by $(n) = (3/8vr)[l - (n-e) 2 ] Q. Note that equation ([[]) has the Lindblad form and, 
thus, Tr[p(t)] = J dzp(z,z,t) — 1. Beside this, for linearly polarized light, P(u) appears to be almost independent of 
u and for most practical purposes one can set P(u) = 1/2. 

To solve Eq. (|l|), we use a Monte-Carlo method of Ref. According to this method one finds the dynamics of 
an arbitrary observable (A(t)) — Tr[Ap(t)] as an average 



(A( t )) = m)\A\m) (4) 

over different solutions of the stochastic Schrodinger equation 

d^ = (~^Hdt ~ \L\L u At + ^L u d^j i> . (5) 

In the latter equation 7 = j(£l/5) 2 is the spontaneous emission rate and d£ is a Wiener process with d£ = and 
d£ 2 = dt (i.e., £(t) is 5-correlated white noise). It is worth to note that the choice of the stochastic Schrodinger 
equation is not unique and one can construct many different stochastic equations corresponding to the same master 
equation. Most often the nonlinear versions, which preserve the norm of the wave function ijj(t), are employed [^5| . 
Here, following Ref. pi , we use a linear stochastic equation which is preferable from the numerical point of view. 
Then the norm is conserved only after averaging over different realizations of the random process £(t). This fact was 
used to control the statistical convergence. 



III. TIGHT-BINDING APPROXIMATION 

First we shall analyze the problem in the tight-binding approximation. In this approximation the Hamiltonian (|^) 
is substituted by the tri-diagonal matrix 

A 

HLm = — — (Sl,m+1 + <>l,m-l) + dFlSi^m , (6) 

and the recoil operator (^) by the diagonal matrix 

Li, m {u) = exp(iirul)8i, m . (7) 

Here d = 7r/fcz, is the period of the optical potential and the index I refers to the localized Wannier function |/) 
associated with the Z-th well of the periodic potential. In the absence of spontaneous emission (7 = 0) and static 
force (F — 0) the eigenfunctions of the system are Bloch waves \ip K ) = J^i exp(id,Kl)\l) corresponding to the energy 
e(n) — — Acos(<i/i). This dispersion relation is assumed to approximate the dispersion relation of the atom in the 
ground Bloch band. In the presence of a static force, the quasimomcntum k of the wave function \ip K ) evolves 
according to the classical equation dn/dt = F/h. In terms of the atomic velocity v — hT 1 de(n) / 'Ok and coordinate z 
this corresponds to a periodic oscillation of the wave packet with the Bloch frequency lob = dF/h: 

(v(t)) = ^ sm{LO B t) , (*(*)) = ^ cos(w B <) . (8) 

(The coordinate and velocity operators are obviously given by the matrices Zi jTn = dl5i, m and w/. m = (dA/2h)(iSi :m+ i — 
i8i >m -\).) Note that within the tight-binding approximation the atoms always oscillate according to a cosine law. In 
reality, however, the oscillations are asymmetric. Besides this, there is a decay of the oscillations due to the interband 
Landau-Zener tunneling |l^, 17, [l8| |l9| - a phenomenon completely ignored by the tight-binding (and more general 



single-band) models. With these remarks reserved, we proceed with the analysis of the effect of spontaneous emission. 

As follows from the explicit form of Eq. ([5]), the recoil operator randomly changes the atomic quasimomentum. 
Thus any narrow distribution of the quasimomentum (which is usually considered as a prerequisite for observing 
Bloch oscillations) will be smeared over the entire Brillouin zone and the oscillations should decay. The following 
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simple "classical" model helps to understand the details of the decay process. Let us consider an ensemble of classical 
particles with the Hamiltonian 



H = - A cos(dp/h) + Fz 
affected additionally by noise (random recoil kicks) so that the conditional probability is 

[d(p-p' - Ft)/h] 2 ' 



W(p,p',t) = (27T7fr 1/2 exp 
Then the mean velocity (v(t)) = (dH/dp) decays as 



4 7 i 



(v(t)) = Vq exp(— jt) sm^LOgt) , vq = dA/h 



(9) 



(10) 



(11) 



and the mean squared velocity relaxes to v 2 t = Vq/2. 

The solid line in Fig. [I] shows the dynamics of (v(t)) calculated on the basis of equations (|5|) - (Q), where as an 
initial condition we choose p(0) = \ip(0))(ip(0)\ with ijj(0) ~ exp(— Z 2 /100) |Z) and set A = 1, h = 1, and d — ir for 
simplicity. It is seen that the solid line in Fig. [I] closely follows the theoretical prediction (11) (dashed line), i.e. the 
decay rate of Bloch oscillations is given by the rate 7 of spontaneous emission. Simultaneously with the decay of the 
mean velocity, the mean squared velocity relaxes to its equilibrium value v 2 t = v 2 /2, which corresponds to a uniform 
distribution of the atoms over the Brillouin zone (so-called recoil heating, see inset in Fig. |l|). 

Since Bloch oscillations decay after a transient time t ~ I/7 it might be naively thought that in the stationary 
regime £ ^ I/7 there is no difference between the cases f ^ and F = 0. Although this is true for the momentum 
distribution, the difference appears when we analyze the atomic dynamics in coordinate space. Indeed, according to 
the "classical model" (^ the velocity correlation function i?t(r) = (v(t + T)v(t)) does not depend on t in the stationary 
regime and obeys 



R(r) 



. exp(— 71-) cos(wbt) 



Substituting Eq. (|12|) into the equation for the mean squared displacement, 



r t \ 2 



(^(*)> = (( / v(t')dt') )^ 2t ] o ^) dT ' 



(12) 



(13) 



we obtain 



(z 2 {t)) - Dt , D = 2v 



st 2 



l 2 



(14) 



Thus the static force suppresses diffusion caused by the recoil heating. As an example, Fig. shows numerical results 
(obtained on the basis of the tight-binding model) for the mean squared displacement (A?(i)) = {z 2 {t)) — {z(t)) 2 
for some values of F (solid lines) together with the diffusion law (Az 2 (t)) = Dt (dashed lines) where the diffusion 
coefficient D is obtained from (14). Concluding this section we stress that Eq. ( |l4"| ) refers only to the stationary regime 
t 3> I/7. In the opposite limit the atomic motion is essentially oscillatory (large static force) or ballistic (weak force) 
and (Az 2 {t)) ~t 2 . 



IV. DECOHERENCE BY RECOIL HEATING 



The results reported in the previous section can be also viewed as a decoherence process. Indeed, in the single- 
band approximation the dynamics of the system is characterized by the density matrix p(t) = ^ n Pn,m{t)\m){n\. 
Substituting Eqs. @-(0) into the master equation (||) and taking into account that P{u) « 1/2 we obtain 

' ~ y^(P»±i,m - Pn,m±i) - iu B (n - m)p n>m - 7(1 - S n . m )p 7hm . (15) 

The last term in this equation causes the decay of off-diagonal matrix elements. Thus the density matrix of the system 
tends to be diagonal in the basis of the localized Wannier states. As an illustration to this statement Fig. || shows the 
absolute values of p n ,m{t) as a gray-scaled map for the parameters of Fig. [l] and t = and t = 5Tb- It is seen in the 
figure that at t > I/7 only the diagonal and nearest to diagonal matrix elements p n , n +i = P* n +\ n have non-negligible 
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values. This property of the " stationary" density matrix allows us to obtain Eq. (|14j) of the previous section without 
appealing to the classical model (P)-(|i~0|). The derivation is as follows. 



First, we estimate the off-diagonal matrix elements p n +i. n - The formal solution for these elements reads 
p n +i,n(t) = ^ e^p{-iu B t - *ft) J dt (p n+1>n +i - p n>n ) exp(ioj B t + jt) 



(16) 



(here we neglect all off-diagonal elements except the ones nearest to the main diagonal). Because in the stationary 
regime the characteristic rate of change of the diagonal matrix elements is much smaller than Wg, we obtain 

Pn+iAt) ~ ttt : ■ (17 

2n iujb + 7 

Next, substituting this estimate into the equation for the diagonal matrix elements yields the rate equation 



dp n , n {t) _ f A\ 2 2 7 



Jo {Pn+l.n+l - 2/9 n.n T pn— l,n— 1 ) • (18) 



dt \2h) uj 2 b +j 2 

Finally, approximating (p n+ i >n+ \ — 2p nj „ + p„_i jn _i)/d 2 by d 2 p/dz 2 we end up with the diffusion equation 

dp _ Dd 2 p 
dt ~ JIh 2 ' 



(19) 



with the same diffusion coefficient D as given in Eq. (jl4|). 

To summarize, the recoil heating reduces the coherence length of the atomic wave function to that of the localization 
length of the Wannier function which, in turn, is of the order of the optical lattice period ir/kL- Obviously, under this 
condition a direct observation of Bloch oscillations is impossible. However, they show up indirectly as a "correction" 
to the diffusion coefficient D. 



V. BEYOND THE SINGLE-BAND RESULTS 



We proceed with the analysis of Bloch oscillation on the basis of Eqs. (jl]) - (|J), i.e. beyond the tight-binding 
(single-band) approximation. It is convenient to use scaled variables, where the length is measured in units of the 
laser wave-length (z — > k^z), the energy in units of the recoil energy Er = h 2 k\/2M , and the time is scaled on the 
basis of the recoil frequency (t — > E^t/h). Then the Hamiltonian (g) takes the form 

H =p 2 -Ucos 2 (z) + Fz , (20) 

where p = — id/dz, U — (hil 2 /S)/ Er, and F is the scaled static force. In what follows we restrict ourselves by 
considering two different values for amplitude of the optical potential: [7 = 1, which we shall refer to as case (a), 
and case (b), where U = 4. The band spectrum of the system in these two cases is depicted in Fig. ^. It is seen in 
the figure that for U = 1 there is an essential deviation of the ground-Bloch-band dispersion relation from the cosine 
law, while for U = 4 the relation €q{k) — — Acos(7r«;) holds with good accuracy. The more important difference of 
two spectra, however, is the size of the energy gap between the ground and "first excited" bands. Indeed, according 
to the Landau-Zener theory this gap defines the probability of interband tunneling, which decreases exponentially 
as a function of squared energy gap |lq ]. Thus, in the case (b) the rate of Landau-Zener tunneling is negligible in 
comparison with case (a). 

Figure || shows the dynamics of an atomic wave packet for 7 = 0, U — 1, and F = 0.025. As expected the 
atoms perform a periodic oscillations with Bloch period Tb = 2/F. It is also seen in the figure that during each 
Bloch cycle (when the wavepacket is reflected at its leftmost position) a fraction of probability is "emitted" in the 
negative direction. In terms of the Bloch-band spectrum this effect corresponds to Landau-Zener tunneling between 
the ground and upper bands, discussed above. Experimentally, this phenomenon was observed in Ref. || and its 
complete theoretical description is given in Ref. |l(J by using the formalism of resonance (metastable) Wannier-Stark 
states. Because of tunneling the probability P(t) to find the atom in any finite interval (larger than the amplitude 
of Bloch oscillations) exponentially decreases with time. In what follows, we quantify this process by the increment 
v. which we refer to as the depletion constant [^of . In numerical simulation the depletion constant was found by 
approximating the function 



/207T 
dz|^(2,i)| 2 
-20ir 



(21) 
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by an exponential function. This gives v = 2.5 • 10 4 and v < 10 6 for the cases (a) and (b), respectively. For the 
sake of future reference the solid lines in Fig. ^| show the dynamics of the mean atomic velocity (v(t)), calculated as 

(«(*)> = / dzr(z,t) (-M-j^) VXM) , (22) 



207T 



and the wave packet dispersion (Az 2 (t)) = {z 2 {t)) — (z(t)) 2 . To compensate the decrease of probability, both (v(t)) 
and (Az 2 (t)) are normalized by dividing by P(t). The asymmetry of the velocity oscillations obviously reflects the 
deviation of the actual dispersion relation eo(ft) from the cosine dispersion relation. The dashed lines in the figure 
correspond to the tight-binding approximation of Sec. 3 where, to take into account the asymmetry of the oscillations, 
the tri-diagonal matrix (^) is substituted by few-diagonal matrix with off-diagonal matrix elements given by the 
coefficients of Fourier transform of eo(n). A reasonable coincidence is noticed. 

We proceed with the case of nonzero rate of the spontaneous emission. First we study the dependence of the 
depletion constant v on 7. These studies are summarized in Table |[ It is seen from the table that the spontaneous 
emission strongly enhances the "tunneling" decay of the system. This can be understood by noting that the recoil 
operator may "kick out" the atom from the ground Bloch band. (More formally, given rp(z) belonging to the subspace 
of Hilbert space spanned by Bloch waves with zero band index, L u ip(z) generally does not belong to this subspace.) 
Because the rate of Landau-Zener tunneling increases with the band index, this causes a faster decay of P(t). 

Since the depletion of the ground band introduces an additional decay mechanism, it might be naively expected 
that the decay of Bloch oscillations should be faster than it is predicted by the tight-binding model. However, the real 
situation appears to be inverse. This is illustrated in Fig. ^, where the dashed lines are predictions of the tight-binding 
model and the solid lines show actual behaviour of the system. (The system parameters are the same as in Fig. || but 
7 = 0.01.) The explanation for this effect is given below. As follows from the explicit form of the recoil operator (||), 
it multiplies the Bloch wave by a plane wave with a wave vector varying from zero to that of the size of the Brillouin 
zone. Within the single band approximation this always results in a change of the quasimomentum k of Bloch wave by 
the same value. In the reality, however, only small values of the wave- vector ("weak kicks") cause a simple change of 
the quasimomentum. The large values ("strong kicks"), as explained above, may remove the system from the ground 
band. As a result, the redistribution of the quasimomentum k over the entire Brillouin zone is slowed down and the 
oscillations of the mean velocity (normalized by survival probability P(t)) decays slower. 

The results of a numerical simulation in the case U = 4 are shown in Fig. |[ Peak-like behaviour of the dispersion 
at short times time origin is a short-time transient phenomenon depending on the initial condition. (In principle, 
this peak can be removed by an appropriate choice of the initial wave packet.) Ignoring this transient phenomenon, 
the overall dynamics of the system is now closer to that predicted by the tight-binding model. This is actually not 
surprising, because the tight-binding approximation is more reliable for larger values of the energy gap. 



VI. SUMMARY 



We analyzed Bloch oscillations of neutral atoms in the presence of spontaneous emission. It is shown that random 
recoil kicks (coming from spontaneously emitted photons) cause a decay of Bloch oscillations through the dephasing 
of Bloch waves. The decay of the mean atomic velocity is accompanied by a diffusive spreading of the atoms in 
configuration space. We studied both of these processes by using the tight-binding approximation and by direct 
numerical simulation of the wave packet dynamics. 

Within the tight-binding model, the increment of Bloch oscillations decay is proven to coincide with the rate of 
spontaneous emission 7 which, in turn, is given by the product of the excitation probability of the upper level and 
its natural width. The diffusion coefficient depends additionally on the value of the Bloch frequency lob (which is 
proportional to magnitude of the static force) and is smaller for larger values of the Bloch frequency. Thus, the static 
force suppresses the spreading of the atoms caused by the recoil heating. 

Direct numerical simulation of the system dynamics confirms qualitatively these results of the tight-binding model. 
However, because of the failure of the single-band approximation, the real situation appears to be more complicated. 
It is shown that the recoil kicks depletes the ground Bloch band with a rate depending on the (scaled) spontaneous 
emission rate 7 and the size of the energy gap following the ground band. Surprisingly, the depletion of the Bloch 
bands affects the decay of Bloch oscillations in a counter-intuitive way: it slows down the decay of the oscillations. 
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TABLE I: Depletion constant v for different rate of spontaneous emission 7. The value of the static force is F = 0.025. 



7 


0.0 


0.001 


0.01 


0.05 


U = 1 


2.5 ■ 10 -4 


6.0 ■ 10" 4 


3.0- 10" 3 


1.0 • 10 -2 


[7 = 4 


< 10~ 6 


1.5 • 10" 4 


1.5 • 10~ 3 


0.9- 10" 2 




FIG. 2: Tight-binding model. Dispersion of the wave packet as a function of time for 7 = 0.05. The slopes of the dashed lines 
are given by values of the diffusion coefficient (H). 
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FIG. 3: Absolute values of the density matrix elements p n , m {t) as a gray-scaled map for t = (a) and t = 5Tb (b). The 
parameters are the same as in Fig. 1. 





FIG. 6: Undamped Bloch oscillations. The mean atomic velocity (a) and the dispersion of the wave packet (b) are shown as the 
functions of time. Parameters are the same as in Fig. 4. The dashed lines show the results of the single-band approximation. 
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